The logistic equation and a linear stochastic resonance 

P. F. GORAt 



M. Smoluchowski Institute of Physics and Complex Systems Research Center, 
Jagellonian University, Reymonta 4, 30-059 Krakow, Poland 



We show analytically that a linear transmitter with correlated Gaussian 
white noises displays a stochastic resonance. We discuss the relation of this 
problem to a generalized noisy logistic equation. 



a > 0, is one of the most frequently used, and undoubtedly most successful, 
models of population dynamics. It is now one of the classical examples of 
self-organization in many natural and artificial systems pp. It is virtually 
impossible to list its all applications — for example, the Web search returns 
more than 60 000 links to pages and on-line publications containing the key 
words "logistic equation". 

A natural generalization of the model is one in which the determin- 
istic growth rate, a, is perturbed by a stochastic process: 



Here ^(t) is the stochastic term. Such equation has been first discussed 
by Leung [2] and later by many other authors. In particular, authors of 
Refs. |H1 13) by Borel summing a formal series for the expectation value 
{x{t)), obtained expressions for the nonlinear relaxation time. The noise 
term ^(t) was represented by the Gaussian white noise (GWN) in Refs. [2J|3] 
and by a Gaussian colored noise in Ref. j4j. Surprisingly, closed and math- 
ematically exact expression for the expectation value and the variance of 
x(t) are still lacking; we have not solved this problem, but we will present 
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1. Introduction 



The logistic equation 




(1) 



X = (a + p^{t))x{l — x) 



(2) 
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a heuristics, based on a simple transformation of Eq. ^ that, as our numer- 
ical simulations indicate, can be useful in predicting the behavior of systems 
described by the noisy logistic equation. 

Surprisingly, the problem of the noisy logistic equation is closely related 
to that of the linear stochastic resonance (LSR). The stochastic resonance 
(SR) is an example of the constructive role of a noise, where the noise 
and a dynamical system act together to reinforce a periodic signal 0; see 
also Ref. jSj for a review. We should mention that an aperiodic stochastic 
resonance has also been discovered [7]. The SR now seems to be ubiquitous 
and has been claimed "an inherent property of rate-modulated series of 
events" |H|- There is, however, a long-standing debate whether the SR 
can be at all present in linear systems. This point is important because 
a linear kinetics (relaxation) is the final stage of most dissipative processes. 
It would be particularly important to establish the existence of a LSR in 
systems driven by a GWN because such a noise represents the standard 
equilibrium fluctuations of macroscopic physical systems. The LSR has been 
first reported by Fulihski in Ref. 9 an later discussed in Refs. |lU[llll[T^ll3j . 
However, it was soon pointed out that the LSR vanished in the GWN limit 
jlUl I12j or after averaging over the initial phase of the signal jllj . or was 
restricted to transient times, even though other constructive effects of the 
noise persisted in the asymptotic regime as well Short-lived or phase- 
dependent phenomena cannot be deemed "unphysical" — on the contrary, 
they may be very important in many situations, for example in enzymatic 
reactions in living cells — but there is a general consensus that a fully- 
fledged stochastic resonance should be characterized by a clear maximum 
of the signal-to-noise ratio (SNR) which survives the phase averaging and 
persists for long times, indicating the presence of an optimal noise level. 

The first evidence for a LSR that displays a maximum of the SNR was 
given only recently in Ref. [^. However, we have shown jl5j that this 
particular case of the SR was of a non-dynamical nature as the output of 
the dynamical system discussed in Ref. 14 merely reproduced the spectral 
properties of the input signal, and moreover, the resonant properties of the 
latter were only an artifact of the chosen parameterization of the noises. 
In the present paper we show analytically that a simple linear system with 
correlated multiplicative and additive GWNs and a driving that consists of 
a constant term and a periodic signal displays a maximum of the SNR that 
lacks all the above-mentioned deficiencies: The maximum is present in the 
asymptotic (long time) regime and not only for transient times, it survives 
averaging over the initial phase of the signal, and finally, it is a result of 
a collective action of the system, the noises and the external driving and not 
merely a reproduction of very particular properties of the input sequence. 
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2. The logistic equation 

Assume that the noise term in the perturbed logistic equation ^ is 
a GWN with (^(t)) = 0, {i{t)^{t')) = 5{t - t'). The equation Q can easily 
be solved for every specific realization of the noise. The formal solution to 
Eq. then reads 



x{t) 



1+ i^exp 



-at-p J S,{t')dt' 





(3) 



In most practical situations, one is interested not in specific realizations of 
a stochastic process, but rather in moments of that process, averaged over 
an ensemble. One is tempted, then, to calculate the moments directly from 
Eq. 0, expanding it in a power series: 



m) = 




where (• • • ) stands for an average over the realizations of the noise £,{t). We 
have changed the order of the summation and taking the average in Eq. Q 
and used the well known fact that for a GWN 




We can, however, see that the final series in Eq. @ is clearly divergent and 
even though a series similar to it has been Borel summed and some useful 
information has been extracted from it [3 13 ) conclude that the formal 
solution given by Eq. © offers little insight into the properties of Eq. © . 
Observe that a formal substitution 



1 

y = - 

X 



(6) 
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leads immediately to 



y{t) = IH exp 

Xq 



t 



at-p S,{t')dt' 







(7) 



and if the process y{t) is convergent (see below), after a sufficiently long time 
almost all realization of the process do not differ from unity and the variance 
of y{t) asymptotically vanishes. This observation does not constitute a proof 
that analogous statements can be made about the logistic process x{t), 
too, but provides a useful heuristics about the long-time properties of the 
logistic process. The substitution ® gives also a possibility of a further 
generalization of the noisy logistic equation that will be pursued in the next 
Section. 



3. The linear transmitter 

With the substitution © Eq. (j2j changes into 

y = -{a + pm)y + {a + pm)- (8) 

Note that this substitution does not affect the stable stationary point x = 1 
of Eq. and its basin of attraction, but it breaks for the unstable station- 
ary point X = of Eq. (^. Therefore, the vicinity of the unstable stationary 
point may not be represented correctly by the transformation 

The equation (jH)) closely resembles the equation of motion of a linear 
transmitter with multiplicative and additive noises that has been the subject 
of our recent research I17j ; the principal difference is that the constant 
additive driving was lacking in models discussed there. The equation 
can easily be generalized 



y = -{a + pU{t))y + {b + qUt)), (9) 

where both noises Ca are GWNs, but they are correlated: {S,m{t)Ca(t')) = 
c6{t — t'). We assume that the noises ?a jointly form a two-dimensional 
Gaussian process (this assumption simplifies the discussion but is not a cru- 
cial one PT|): 



(lOa) 
(10b) 
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where ^, r/ are uncorrelated (independent) GWNs and c G [—1, 1]. We finally 
arrive at 



-(a + pi{t))y + h + qci{t) + qVT^ r]{t) . 



3.1. Moments of the process y{t) 
A formal solution to Eq. ((TT|l reads (y(0) = 0) 



e exp 



-P J i{t")dt" 
t' 



b + qc^it') + qVl-c'^vit')) dt' . 

(12) 

Unlike the solution ©, Eq. (O can be used directly to calculate the 
moments: 



{y{t)) = b / e-'^(*-*') 



exp 



-P I ^dt" 



dt' 



+ qVl - c2 J e-°-^'~''Uri{t')ex.^ -p j i{t")dt 



dt' 



+ qc J e 




-a{t-t') 



i{t') exp 



-p j i{t")dt" 



dt' . 



(13) 



The first of the expectation values on the right hand side of Eq. (|13j) is of 
the form (jSJ. The second vanishes identically because the noises r] and ^ 
are uncorrelated, and we have calculated the last one in Ref. |17j : 



C(t')exp -pj^{t")dt" 
t' 

After collecting all terms we get 



^P exp Q/(t-t')) . 



(14) 



(yit)) 



1 9 

a — ^p'^ 



1 - e" 



h 1- 



^ - (, „-(a-ip2)A ^ _ - jcpq 

t-»oo a — 72P 



(15) 



The limit in Eq. (|15|) exists if a — > 0; otherwise the process is 
divergent. If we want to keep the interpretation that rr = 1/y is a measure 
of a population, the limit in Eq. ()15|) makes a physical sense if 6 — ^cpq ^ 0. 
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We now calculate the variance of the process y{t): 



t t I L L 

{y\t)) = JdhJ dt2e-'^(2t-ti-t,) / _p I ^{t')dt'-p J i{t')dt' 

^ L tl t2 . 



b + qcCih) + qy^l^7]{ti)) (b + gcCfe) + g\/T^r/(t2) 



(16) 

As we can see, one needs to know one more expectation value that we have 
previously calculated in Ref. 



exp 



fimt')dt' 



Lo 



[Sih - h) + f{ti)f{t2)] exp 



T 



[m]'dt' 



(17) 



It is now easy to verify that the variance (|16|) exists if a — > 0. The full 
expression for the variance is rather long, and therefore we present the final 
expression in the asymptotic regime only: 



D = {y\t)) - {yit)Y 



t~*oo 



Ab'^p'^ - Sabcpq + {4a^ - 4a(l - c^)p'^ + (1 - c^)p'^) 
2(a-p2)(p2 _2a)2 



(18) 



Let us consider the limiting cases. First, for the uncorrelated noises we 
obtain 



D 



b'^p^ + (a - \p^Yq^ 



c = 0. 



(19) 



2(a-p2)(a_ ip2)2 ' 

As we can see, for the uncorrelated noises D increases monotonically both as 
a function of the amplitude of the multiplicative noise, p, and as a function of 
the amplitude of the additive noise, q. Next consider the case of maximally 
correlated (or maximally anticorrelated) noises: 



D 



{bp =F aqf 



2(a-p2)(a- ip2)2 ' 



c = ±1. 



(20) 
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The variance (|2()|) can vanish for certain values of parameters: D = if 
c = ±1 and bp^ aq = 0. For such parameters, almost all realizations of the 
process y{t) asymptotically remain equal their expectation value yoo despite 
the fact that the noises do not cease and constantly act on the system. 
Note that in this case yoo = ^q/p the negative value corresponds to an 
"unphysical" realization of the process (see the discussion below Eq. 

In general, the variance D exhibits, as a function of a minimum for 
all |c| > and appropriate values of other parameters; some of these sets 
of parameters may correspond to "unphysical" realizations. D is, by con- 
struction, nonnegative and equals zero only at its minimum for c = ±1. 

The fact that the linear transmitter described by Eq. admits, for 
the maximally (anti) correlated noises, realizations with a vanishing variance 
can be shown quite elementarily: Let c = ±1 and bp^aq = 0. Then Eq. pip 
is equivalent to 



y = -{a + pm)y + {b±qm) = -{a + Pm){yT 



p 



y{t) = + (2/(0)t - )exp 
p V P, 



-at-p j (,{t')dt' 



(21) 



(22) 



If the process is convergent, almost all realizations asymptotically approach 
Voo = M/p- 



3.2. Application to the logistic equation 

According to what we have said, the linearized logistic equation (jSJ 
clearly leads to a solution with asymptotically vanishing D, or with (^y'^{t)^—f 

{y{t))'^ — > 1. While in general there is no simple relation between the vari- 
ances D = (y^(t)) — {y{t))'^ and = (x^(i)) — {x{t))'^, we expect that if 
D = and yoo ^ 0, Dx vanishes, too. Observe that for sufficiently large 
times, with (a — |p^)~^ being the characteristic time scale, almost all real- 
izations of y(t) reach a constant value, and so does the inverse l/y{t) = x{t). 
We thus conclude that almost all realizations the process described by the 
logistic equation with a noisy growth rate © reach a constant value, equal 
to the equilibrium value of the deterministic logistic equation Heuris- 
tically this means that once the process reaches a value x{t) ~ 1, the right 
hand side of Eq. ^ remains nearly equal to zero for all future times even 
though the growth rate fluctuates. We have confirmed these results by direct 
numerical simulations. 
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Fig. 1. Average distributions of values of x{t) described by Eq. H23(l as a function 
of the additive noise level, q. Parameters are: a — 2, b — 1, p = 0.125. The 
correlations are c — I (upper panel) and c — 0.5 (lower panel). 



The equation for the Hnear transmitter corresponds to a further 
generalization of the logistic equation 

x = -{a + p^{t))x + {b + qc^{t) + qVl-c^vit))x'^ , (23) 

where not only the growth rate, but also the limiting population level fluc- 
tuate, and these fluctuations are correlated. The results of the previous 
subsection suggest that in the presence of correlations, for certain values of 
the parameters the variance of x(t) should first shrink and then grow as the 
level of the additive noise increases. We have performed numerical simula- 
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tions to confirm this. We have solved the nonhnear equation (|2I-{|) using the 
Heun scheme |18| with a time step h = 1/256. The GWNs have been gener- 
ated by the Marsaglia algorithm and the famous Mersenne Twister |2U| 
has been used as the underlying uniform generator. We have let the system 
to equilibrate for 2^^ time steps, collected the results for the next 2^^ steps, 
calculated the distribution of the values of x{t) that have appeared during 
the simulation, and averaged the results over 128 realizations of the process. 

Selected results are presented in Fig. ^ In case of maximal correlations 
(c = 1, upper panel) the distribution of the values of x{t) first gets narrower 
as q increases, becomes practically (5-shaped for q = hp /a, and then widens 
as q increases past the resonance. A similar, but much weaker, effect is 
observed in case of partial correlations (c = 0.5, lower panel). Note that for 
large values of q the distributions become skewed, with a marked preference 
for values above the expectation value. 



So far we have discussed the systems without any external signal acting 
on them. Now consider an additively coupled periodic signal: 



ij = -{a+pi{t))y + b + qci{t) + q^J I - r]{t) + Acos{nt + (/>) , (24) 



where cp is the initial phase of the signal, A is the amplitude, 0, is the 
frequency, and all the other parameters are as above. We will, in addition 
to taking the usual average over realization of the noises, average over the 
initial phase, as otherwise the correlation function would not correspond to 
a stationary series pi] : 



The calculations leading to the final formula for the correlation function 
are straightforward but tedious. The expectation values ©, (EJ), and 
H17|) are the most important ingredients. We obtain 



4. A linear stochastic resonance 




(25) 







{{yit)y{t + T)))-{{y{t))f- 




+ D e-("-^p')^ 



2[0 



(26) 
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0.10 0.12 0.14 



Fig. 2. Signal-to-noise ratio (SNR) for the linear transmitter (|24ll with various 
correlations between the noises: c = 1 (solid line), c = 1/2 (broken line), c = 1/4 
(crosses). Amplitude of the signal A = 1, frequency fl — 2tt. Other parameters as 
on Fig.^ Curves with c' — — c and 6' = —6 are identical to those presented above, 
cf. Eq. (UHl). 

where D is given by Eq. H18|) . We can now calculate the power spectrum of 
the process y{t) as the Fourier transform of the above correlation function 
(Wiener-Khinchin theorem) and then the signal-to-noise ratio (SNR) as the 
ratio of the power associated with the signal, which corresponds to the 
Fourier transform of the first term on the right hand side of Eq. ()26() , to the 
power of the noisy background at the frequency of the signal: 



SNR = 10 log 



10 



signal 



P.. 



(27) 



By convenience, the SNR is usually measured in dB. 
For c = 1 this procedure yields 



SNR = 10 log 



2A'^{a-p'^){a 



10 



ip2)2 + 2 [(C 



aqf 



(28) 
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As we can see, the SNR has, as a function of g, a clear maximum for bp—aq = 
0. The presence of this maximum provides an unequivocal evidence for the 
stochastic resonance. It is easy to verify that for all [c[ > and appropriate 
values of other parameters the SNR has a maximum as a function of the 
additive noise strength, but for small values of [c[ this maximum is not very 
much pronounced, cf. Fig. [2 Note that also away from the maximum, larger 
values of |c| correspond to larger values of the SNR. 

We have confirmed these results by direct numerical simulations. Nu- 
merical results also indicate that there is a stochastic resonance in case of 
a multiplicatively coupled signal. Analytical results are difficult to obtain 
in this case due to the presence of transcendental function in the expression 
for the correlation function, cf. Refs. |16[ I17j . 

5. Discussion 

A clear and conclusive evidence that the linear system (|24|) displays 
a fully-fiedged stochastic resonance when driven by Gaussian white noises 
is the principal result of this paper. The linear stochastic resonance re- 
ported here is characterized by a clear maximum of the SNR, persists for 
asymptotically long times and survives averaging over the initial phase of 
the signal. This result closes, we believe, the long debate whether the linear 
stochastic resonance is at all possible. It is important to realize that there 
are two factors that are needed to produce the SR in the system (|24|) : (i) 
the multiplicative and additive noises have to be correlated, and (ii) there 
should be a constant driving term, b ^ 0, present. Berdichevsky and Git- 
terman in Ref. [12^ and the present author in Ref. |17| have also considered 
correlated multiplicative and additive noises, but without the constant driv- 
ing, no stochastic resonance has been present, at least for the GWN case. 
It is worth mentioning that the authors of Ref. [T3| also discussed a system 
with two correlated Gaussian noises, but in that case one of the noises was 
coupled multiplicatively to the signal. 

Note that the system ()24() is linear in the sense that its equation of 
motion is linear with time-dependent coefficients, and thus we call the phe- 
nomenon discussed above a "linear stochastic resonance." There is, how- 
ever, a point in observation that since noise is meant to represent many 
unobserved and unaccounted for degrees of freedom — instead of consider- 
ing impossibly complicated microscopic motions, we mimic their effect by 
a reasonably simple stochastic process — a multiplicative (nonlinear) cou- 
pling between the stochastic process and the observed degrees of freedom 
means a "hidden" nonlinearity. 

We have also shown that the problem of the linear transmitter is closely 
related to the generalized noisy logistic equation. The analytical results for 
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the transmitter provide some heuristics that might be helpful in predicting 
the behavior and properties of the noisy logistic process. However, the prob- 
lem of finding closed and mathematically exact formulas for the moments 
of the latter remains open. 

I am grateful to Dr. Ryszard Zygadlo for helpful discussions and to Prof. 
Andrzej Fulihski for his constant encouragement. 
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